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Detection thresholds in polarized intensity and polarization bias correction are investigated 
for surveys where the polarization information is obtained from RM synthesis. Considering 
unresolved sources with a single rotation measure, a detection threshold of &<Jqu applied to the 
Faraday spectrum will retrieve the RM with a false detection rate less than 10~ 4 , but polarized 
intensity is more strongly biased than Ricean statistics suggest. For a detection threshold of 
Sctqu, the false detection rate increases to ~ 4%, depending also on A 2 coverage and the extent 
of the Faraday spectrum. Non-Gaussian noise in Stokes Q and U due to imperfect imaging 
and calibration can be represented by a distribution that is the sum of a Gaussian and an 
exponential. The non-Gaussian wings of the noise distribution increase the false detection rate 
in polarized intensity by orders of magnitude. Monte-Carlo simulations assuming non-Gaussian 
noise in Q and U, give false detection rates at 8aQu similar to Ricean false detection rates at 
4.9aQij. 
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1 Introduction 

Linear polarization of radio sources contains informa- 
tion on magnetic fields in these sources, and Faraday 
rotation of the plane of polarization provides informa- 
tion on the direction and magnitude of the magnetic 
field along the line of sight. As such, observations of 
linear polarization of radio sources provide the most 
widely applicable probe of cosmic magnetic fields on 
scales from galaxies to clusters of galaxies. Finding 
polarized sources in survey images and fitting their 
parameters forms the basis of this analysis. 

Sources with detectable polarized emission are read- 
ily identified in images of total intensity. As a first 
approximation, source finding in polarization can be 
reduced to applying a suitable detection threshold to 
the polarized intensity at the location of every radio 
source identified in total intensity. In practice, source 
finding in polarization is more complicated because of 
two reasons. The first is related to resolved sources, 
and the second is related to Faraday rotation of the 
polarized emission. 

Figure [TJ shows a radio source from the National 
Radio Astronomy Observatory (NRAO) Very Large 
Array (VLA) Sky Survey (NVSS: ICondon et al. 111998ft 
that is slightly resolved in total intensity (white con- 
tours), with major axis 68" in position angle —72°. 
The polarized emission is shown in grey scales as Stokes 
Q and U images. The source has a component that is 
unresolved at the 45" resolution of the NVSS. The un- 
resolved component has a peak polarized intensity of 
6 mjy (uqu = 0.35 mjy), approximately twice the cat- 
alogued value derived from the polarized intensity at 
the fitted position of the source in total intensity. It is 
not clear which fraction of radio sources has different 
morphology in polarized intensity than in total inten- 



sity. Approximately 8% of NVSS sources brighter than 
10 mjy outside the Galactic plane (\b\ > 30°) has a 
fitted (i.e. before deconvolution) major axis size more 
than 1.5 times the 45" (FWHM) size of the synthe- 
sized beam, and 2% has a fitted major axis more than 
twice the beam size. 

Resolved polarized sources have been treated in 
different ways in the literature. The NVSS catalog 
derived polarized intensity at the location of the fitted 
position in total intensity, and the listed polarized flux 
density (peak times Stokes / solid angle) implicitly as- 
sumes a cons t ant p olari zation angle over th e source. 
iTavlor et ai~l 1120071 ) and lGrant et al. I (|201Ch fitted 2- 
dime nsional Gaussians to sources in polarized inten- 
sity. Sub rahmanvan et al. I (2010) integrated Stokes 
Q and U over the solid angle of the source defined in 
a low-resolution total intensity image, and catalogued 
each polarized source as if it were unresolved. 

The second complication for source finding in lin- 
ear polarization is presented by Faraday rotation, even 
if the source is unresolved. Faraday rotation rotates 
the polarization angle by an amount proportional 
to A 2 for a simple Faraday thin source, and in a more 
complicated manner e.g. if synchrotron emission and 
Faraday rotation both occur in the same volume. Dif- 
ferential Faraday rotation over the observed frequency 
range leads to depolarization, thus introducing a selec- 
tion effect against sou rces with strong Faraday rotation 
(|Stil fc Tavlor 1 120071 . for the NVSS). Multifrequency 
observation s allow Rotation Measure (RM ) Synthesis 
(|Burnlll966l : iBrentiens fc De Bruvnll2005h to solve for 
the unknown Faraday depth, polarized intensity, and 
polarization angle simultaneously. The source is typ- 
ically identified by the maximum value in the Fara- 
day spectrum , e.g. through the RM clean algorithm 
(|Heald Il200"9l h 
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Figure 1: A radio source that is resolved in total 
intensity, with one-sided polarization. Grey scales 
show Stokes Q and U from the NVSS with con- 
tours of total intensity at 3, 30, and 60 mJy. The 
grey scales range from —3 mJy (black) to +3 mJy 
(white). 



Following iBrentiens fc De Bruvr] (|2005h the Fara- 
day depth <f> is denned as 



<p(r)= 0.81 



-Biin e • dr rad m 



(1) 



where Bu is the line of sight magnetic field com- 
ponent, n e is the thermal electron density, dr is an 
infinitesimal path length, with the integral taken from 
the observer to the point x. 

The complex polarized intensity V (A 2 ) = Q+iU is 
the Fourier transform of the Faraday dispersion func- 
tion F ((/)), 



P(A 2 ) 



F^e 2 ^ d<f>. 



(2) 



The Faraday rotation measure RM is defined as 
the slope of a polarization angle versus A 2 plot: 



AM (A) 



where 



• tan 



-i U 



(3) 
(4) 



Once a polarized source has been detected, the 
observed polarized intensity p = y/ Q 2 + U 2 must be 
corrected for polarization bias. Since p is a positive- 
definite quantity, the noise in Stokes Q and U results in 
a positive value for p even if no signal is present. The 
statistics of p with a si gnal pp and n oise oqu is given by 
the Rice distribution ( Rice I fl~945h . Estimators of the 
true polarized flux density po from the observed polar- 
ized flux density p based on this d istribution have been 
discu s sed by a number of authors ((Simmons fc Steward 
119851 : I Vaillancouril 120061 ). The Rice distribution as- 
sumes Gaussian noise in Stokes Q and U. Real surveys 
of the sky may have non-Gaussian tails to the noise 
distribution resulting from imperfect imaging and cal- 
ibration. In this paper we investigate the detection 
statistics in polarized intensity and polarization bias 
correction for polarized intensity determined from RM 
synthesis and in the presence of non-Gaussian noise in 
Stokes Q and U . 



Future radio polarization surveys such as the Galac- 
tic Arecibo L-band Feed Array Continuum Transit Sur- 
vey (GALFACTS; iTavlor fc Salter ll201Ch and the Po- 
larization S ky Survey of the Univ erse's Magnetism 
(POSSUM; iGaensler et al~ll2010r) are wide band, mul- 
tifrequency surveys that require revision of detection 
threshold and polarization bias correction. 



2 RM Synthesis simulations 



Table 1: False detection rates in simulated Fara- 
day spectra. (1) Ratio of true polarized signal to 
the noise in Stokes Q and U. (2) Percentage of sim- 
ulated spectra with | </> — 0o | > 50 radm -2 . No 
false detections were found in the 3.33 x 10 4 8a qu 
simulated spectra. These percentages depend in 
part on the Faraday depth range considered, and 
the width of the RM spread function. 



Po/o-qu 

(1) 


% False Detection 
(2) 


2.0 


73.9 


3.0 


43.9 


4.0 


16.7 


5.0 


3.6 


6.0 


0.43 


7.0 


0.033 


8.0 





RM synthesis was performed on simulated data for 
sources with signal to noise ratio po/ctqu ranging from 
to 15. Each realization contained a source with the 
prescribed polarized signal and random polarization 
angle at the reference frequency, 1400 MHz. Stokes 
Q and U values were calculated assuming a Faraday 
depth of 150 radm -2 for A^an = 1024 frequency 
channels between 1000 MHz and 1400 M Hz. G aus- 
sian noise with standard deviation a — V^chan was 
added to each channel, resulting in a standard devia- 
tion (jqu = 1 for the noise in Stokes Q and U after 
averaging over all channels. RM synthesis was then 
performed on the synthetic spectrum of complex po- 
larization. The effect of spectral index a (S„ ~ 
assuming that the percentage polarization is constant 
across the band, was investigated with separate sim- 
ulations for a — and a = —0.75. For each com- 
bination of Po/oqu and a, 33 300 simulated Faraday 
spectra were analyzed. 

Figure [2] shows two simulated Faraday spectra with 
and without noise. The noiseless spectra show the RM 
spread function with its side lobes. The near side lobes 
of the RM spread function raise the probability of a 
false peak at the wrong Faraday depth, resulting in 
stronger wings in the error function of the Faraday 
depth when the signal to noise ratio is low. As the 
Faraday depth of the source is not known a priori, the 
location of the peak in the Faraday spectrum is cat- 
alogued as the Faraday depth of the source, and the 
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amplitude of the peak as the observed polarized inten- 
sity p m ax- We also extract the polarized intensity p at 
the input Faraday depth, because of the expectation 
that the Rice distribution with noise <tqu = 1 applies 
to p, not necessarily to p max . 
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Figure 2: Faraday spectra of simulated sources 
with signal to noise ratio Po/tqu = 5 (a) and 
10 (b), both with Faraday depth of 150 radm~ 2 . 
The red curves represent the spectrum with noise, 
while the blue curves show the corresponding spec- 
trum without noise. 

Figure [3] shows the distribution of Faraday depths 
derived from the simulations for p/oqu = 2 and 4 
and a = 0. While each Faraday profile contains a 
source with the input 4>o = 150 radm -2 , the uniform 
distribution of Faraday depths at low signal to noise 
ratios represents spurious peaks in excess of the actual 
source. The effect of these false detections is twofold: 
expectation value of the polarized intensity approxi- 
mates a nearly constant value for Pq/<jqu < 3, and 
the error distribution in Faraday depth becomes sig- 
nificantly non-Gaussian. Table Q] lists false detection 
rates as a function of signal to noise ratio. In the range 
3 < Po/oqu < 8 the false detection rate drops to be- 
low 10 -4 , and the error distribution of Faraday depth 
becomes approximately Gaussian. 

Figure U shows the distributions of polarized inten- 
sity derived from the simulations, along with curves 
representing the Rice distribution for Pq/oqu — 2, 3, 
4, and 5. The red histogram shows the polarized in- 
tensity p$ at the actual Faraday depth of the source, 
demonstrating the Ricean statistics for this quantity. 
In reality, we do not know the actual Faraday depth of 
the source, but solve for this by finding p ma x, defined as 
the peak of the Faraday spectrum. The distributions 
of Pmax are shown by the blue histograms in Figure [4] 
At any signal to noise ratio, the distribution of p ma x is 
shifted to higher p with respect to the Rice distribu- 
tion, adding to the polarization bias. This additional 
bias is closely related to the fitting bias for fitting the 



Figure 3: Distribution of Faraday depths, with in- 
put 0o — 150 radm -2 , derived from simulations 
with (a) p/o-qu = 2 and (b) p/<T QU = 4. 



flux d ensity of a source discussed by ICondon et al. I 
(I1998T 1. The magnitude of this bias is similar to the 
well-known polarization bias. 

For po/oqu < 5, the distribution of p max is skewed 
with respect to the Rice distribution. The distribu- 
tion of p m ax can be approximated by the distribution 
of the maximum of TV — 1 independent draws from the 
Rice distribution with no signal, and 1 draw from the 
Rice distribution with the seeded signal po, where TV 
is the ratio of the range of the Faraday spectrum to 
the FWHM width of the RM spread function. The 
sidelobes of the RM spread function increase the false 
detection rate, so this explanation can only be an ap- 
proximation. The difference between the distribution 
of Pmax and the Rice distribution at low signal to noise 
depends on the range of the Faraday spectrum. At 
high signal to noise ratios, the maximum is always as- 
sociated with the source, and the fitting bias is inde- 
pendent of the range of the RM spectrum. 

Table [2] lists expectation values of polarized in- 
tensity for flat spectrum sources and steep spectrum 
sources at a range of signal to noise ratios in polarized 
intensity. Column 6 lists the estimated true polarized 
intensity using p^ = yjv\ a ~ a Qv Though this is not 
measurable in real data, the correspondence between 
column 6 and column 2 reflects the Ricean statistics of 
p</, illustrated in Figure 0] The bias inp max measured 
from real data is approximately twice as large as in 
p</, . The polarization bias correction can be adjusted 
to correct for the additional bias associated with the 
uncertainty in RM. The estimator 

po = ^-2.3(7^ (5) 
for p/oqu > 4 provides accurate estimates of po. 
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Figure 4: Distribution of polarized intensity de- 
rived from simulated Faraday spectra. Signal to 
noise Po/o-qu = 2, 3, 4, 5 for panels (a)-(d) respec- 
tively. The red histograms represent values of p at 
the actual Faraday depth of the source, and the 
black curves show the Rice distribution for the as- 
sumed value of Po/o-qu- The blue histograms rep- 
resent polarized intensity at the peak of the Fara- 
day spectrum. 



This work suggests that a detection threshold of 
Po/pqu > 8 should be applied for the derivation of 
Faraday depth, in order to obtain a well-behaved error 
function of Faraday depth. Polarized intensity can be 
estimated for sources with Po/ctqu > 4, depending on 
the desired level of acceptable false detections. 

In the case where a 7^ 0, [Brent i ens fe De Bruvnl 



(2005) recommend dividing by the total brightness as 
a function of frequency. This may work well for bright 
sources, but not for faint sources, or diffuse polarized 
emission. In our experiments, we find that the Faraday 
depth derived for sources with a = —0.75 was not 
significantly different from sources with a = 0, but the 
polarized intensity after polarization bias correction is 
given by p e f j defined as 



Peff 



J pd\ 2 



(6) 



where the integral is evaluated over the wavelength 
range of the data. The penalty of not dividing by 
total intensity creates a spectral-index dependent bias 
in polarization that is larger than the effects discussed 
previously. 
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Figure 5: Example of one of the 4° x 4° fields in 
(top) total and (bottom) polarized intensity, illus- 
trating the effects of missing fields on the noise. 



3 Simulated Sky Survey 

Sky simulations with sensitivity and angular resolution 
similar to the NVSS over an area of 2.33 sr covering 
478 4° x 4° fields were constructed to test source find- 
ing and stacking algorithms. Figure [5] shows a simu- 
lated image in total intensity and polarized intensity. 
The images were constructed with source density and 
noise level similar to the NVSS polarization images. 
The images were built up by seeding sources at ran- 
dom positions in the image plane. Each source con- 
sists of a VLA snapshot antenna pattern scaled to the 
assumed clean limit plus a two dimensional Gaussian 
representing the restored clean components. The side- 
lobes of the antenna patterns from different sources 
add up incoherently simulating incomplete cleaning. 
The Gaussian noise divided by the sensitivity pattern 
of the NVSS mosaics was added to the Stokes /, Q and 
U images. The noise level in Stokes / is 2 mjy, whilst 
Q and U is 0.4 mjy. Images with missing fields were 
included to simulate survey edges. The distribution of 
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Figure 6: Comparison of the ratio of peak values 
(J and p) retrieved from the image and the input 
values (Iq and po) for sources with p > 20 mJy. 
(a) shows peak fluxes determined from fitting with 
SAD, (b) shows the nearest pixel value. 



polarized f raction of sources f ollowe d the distribution 
derived bv lBeck fe Gaensler I (|2004ft for NVSS source 
brighter than 80 mJy. The images only contain unre- 
solved sources and the resolution of the images is 45" 
with a pixel size of 15". 

The images were searched for sources in total in- 
tensity and polarized intensity with the AIPS (Astro- 
nomical Image Processing System) source finder SAD 
(Search and Destroy), with a detection threshold of 
5 mJy. The recovered sources were matched with the 
input source catalogue. Only sources that matched 
within 60" were considered for further analysis. For 
sources with I > 50 mJy the standard deviations in 
right ascension and declination were 1.47" and 1.46" 
respectively. 

The peak flux in both total and polarized intensity 
is found in two ways. By using SAD to find and fit the 
sources and by extracting the nearest pixel value for / 
and p. The nearest pixel values are considered because 
RM synthesis is done per pixel. Figure [6] compares the 
fitted peak and the nearest pixel values with the input 
catalog for both total intensity and polarized intensity 
for sources with p > 20 mJy. The fitted peak values 
underestimate the input values by a few percent, but 
the solid angle of the source from the fits is slightly 
overestimated so that the integrated flux density is re- 
trieved from the catalogue. The nearest pixel inten- 
sities underestimate the true intensity by up to 15%, 
approximately along the line of constant p/I. The er- 
ror in the fitted position of the source is much smaller 
than a pixel, so the uncertainty in the position of the 
source in total intensity does not introduce a signifi- 
cant error in the estimation of p. RM synthesis on the 
brightest pixel would introduce a systematic error in 
polarized intensity comparable to that shown in Fig- 
ure [6b. Figure [6] suggests that source finding in the 
image plane after RM synthesis is required to obtain 



polarized flux densities with an accuracy better than 
~ 10%. 

4 Non-Gaussian noise in Q/U 

The statistics of polarized intensity are usually de- 
scribed by the Rice distribution that assumes Gaussian 
noise in Q and U. Imperfect imaging and calibration 
result in images that do not achieve the theoretical 
noise levels. The actual noise distribution has strong 
wings above Gaussian noise. 




Stokes Q Flux (mJy) 



Figure 7: Histograms of the noise for Stokes Q im- 
ages for (a) the simulated images (b) the NVSS. 
The fitted curves are (solid) a Gaussian fit and 
(dashed) the summation of a Gaussian and an ex- 
ponential. 

Figure [7] shows the noise distribution for both the 
simulated and the NVSS images. Sources identified in 
total intensity were masked out of the Stokes Q and U 
images leaving pixels that are free of detectable polar- 
ized emission. For determining the noise statistics, ar- 
eas near missing fields were not considered. The solid 
curves in Figure [7] represent Gaussians with standard 
deviation equal to the rms of empty areas in the im- 
ages avoiding sources. The non-Gaussian wings in Q 
and U emerging above the Ioqu level are related to 
the striping visible in Figure [5] In both cases a Gaus- 
sian does not adequately represent the wings of the 
noise distribution. A better solution is the sum of a 
Gaussian and an exponential, 

F{x) = Ae-* 2/2CT2 +Be~ CM , (7) 

with A, B and C determined from fitting. The pa- 
rameters of the fit for the Stokes Q simulated images 
are: A = 0.95341, B = 0.00659, C = 1.08135 and 
a = 0.44327. 

To investigate the impact on false detection rate in 
polarized intensity, Monte-Carlo simulations of polar- 
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ized intensity using the noise distribution from Equa- 
tion [7] and po — were done. First a Q value was 
drawn. In principle Q values were drawn indepen- 
dently from U values. However, if the \Q\ was larger 
than a Igqu threshold, U values were drawn until the 
\U\ was also larger than the 2oqu- This procedure ac- 
knowledges that in real data residual sidelobes in Q 
probably also exist in U . 

Figure [8] shows the false detection rate in polar- 
ized intensity for Gaussian and non- Gaussian noise in 
Stokes Q and U and Table [3] lists false detection rates 
for a range of detection thresholds pu m in polarized in- 
tensity. The details of Figure[S]depend on the dynamic 
range of the Q and U images and will be different for 
every survey. The non- Gaussian wings in Q and U in- 
crease the false detection rate by orders of magnitude. 
In the simulations an 8cr detection threshold yields the 
same false detection rate as a 4.9cr detection threshold 
for Ricean statistics. 

Polarized source finding should apply a detection 
threshold that is derived from the actual noise distri- 
bution in Q and U of a dynamic range limited survey. 
We found no significant effect on the bias correction 
from the non-Gaussian wings. 



Table 3: False detection rates of polarized intensity 
for Gaussian (i.e. Ricean) and non-Gaussian (i.e. 
non-Ricean) noise in Stokes Q and U. 



Piim/o~Qu Ricean Non-Ricean 



3.0 


1.36 x 1Q- 


i 


1.41 


X 


10 


-i 


3.5 


4.43 x 10- 


2 


4.88 


x 


10 


-2 


4.0 


1.13 x 10- 


2 


1.51 


x 


10 


-2 


4.5 


2.23 x 10- 


3 


5.34 


x 


10 


-3 


5.0 


3.35 x 10- 


4 


2.82 


x 


10 


-3 


5.5 


3.85 x 10" 


-5 


2.00 


x 


10 


-3 


6.0 


2.90 x 10- 


-6 


1.54 
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10 


-3 
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3.00 x 10- 
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1.20 


x 


10 


-3 


7.0 
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x 


10 
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x 
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-4 


8.0 
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Figure 8: False detection rates determined from 
Monte-Carlo simulations of the noise distribution 
for po — 0. The dashed curve represents a Gaus- 
sian fit to the noise and the solid curve uses the 
sum of the Gaussian and exponential noise func- 
tion. 



5 Summary &; Conclusions 

The uncertainty in the Faraday depth of the source in- 
troduces a stronger bias in polarized intensity than just 
the well known polarization bias. At low signal to noise 
the false detection rate is greatly enhanced, while at 
higher signal to noise (p > 4(tqu) an effective estimator 
for the true polarized intensity is po = ^Jp 2 — 2.3a q V . 

RM synthesis on the pixel nearest to the fitted po- 
sition of total intensity introduces a systematic error 



that underestimates the polarized intensity by up to 
15%. Source fitting in polarized intensity provides a 
more accurate result, even for the unresolved sources 
considered in this paper. 

Non-Gaussian wings of the noise distribution in 
Stokes Q and U significantly increase the rate of false 
detection in polarized intensity by orders of magni- 
tude. False detection rates at 8<jqu similar to Ricean 
false detection rates at 4.9<jqu- 
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Table 2: Expectation values of polarized intensity and bias correction. (1) Spectral index (S v ~ v a ). (2) 
Input polarized intensity at centre of the band in units of o~qu, and constant percentage polarization across 
the frequency band. (3) Effective polarized intensity defined in Equation 2. (4) Polarized intensity in the 
Faraday profile at the input Faraday depth (150 radm" 2 ). This quantity is not known for real sources. (5) 
Maximum polarized intensity taken over all Faraday depth values. (6) Estimator of Po,eff taking polarized 

intensity from column 4, according to p<f, = Jv\ — a hu- This quantity is not known for real sources. 



Dr of p . c 


// taking 


polarized intensity from column 5, according to po 


Y Pmax 
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Pn/o-Qu 


P0,eff/<?QU 


P<Po/ a QU 
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4.110 


4.246 


4.467 


4.130 


4.202 


-0.75 


5.000 


5.138 


5.244 


5.365 


5.150 


5.147 


-0.75 


6.000 


6.165 


6.255 


6.345 


6.175 


6.161 


-0.75 


7.000 


7.193 


7.272 


7.347 


7.204 


7.188 


-0.75 


8.000 


8.221 


8.288 


8.353 


8.227 


8.215 


-0.75 


9.000 


9.248 


9.305 


9.362 


9.251 


9.238 


-0.75 


10.000 


10.276 


10.332 


10.383 


10.284 


10.272 


-0.75 


15.000 


15.413 


15.444 


15.478 


15.412 


15.403 



